Finite-width QCD sum rules for p and uo mesons* 



B. Kampfer, S. Zschocke 



Abstract 



Forschungszentrum Rossendorf, PF 510119, 01314 Dresden, Germany 

g ; February 8, 2008 

O 

> 
O 

. We present a combined analysis of the in-medium behavior of p and u mesons within the Borel 

QCD sum rule taking into account finite widths. 

(N 1 Introduction 

The experiments at the High-Acceptance Di-Electron Spectrometer HADES pQ are aimed at verifying 
predictions of the behavior of light vector mesons in nuclear matter. Due to the decay channel V —>■ e + e~ 
g \ a measurement of the escaping di-electrons can reveal directly the properties of the parent mesons, 
■ V — p, u, ■ ■ •, since the interaction probability of the e ± with the ambient strongly interacting medium 
+-p> ■ is small. 

In relativistic heavy-ion collisions temperature effects play an important role, which also cause a 
^ ! change of the properties of vector mesons. Indeed, the experiments of the CERES collaboration at 
" the CERN-SPS [I can only be explained by assuming strong medium effects, in particular for the p 
meson (cf. |H E] and further references therein). This seems to be confirmed at higher beam energies, 
^ ■ as delivered by the Relativistic Heavy Ion Collider, since the p meson, as measured via the 7r + 7r~ decay 
channel, suffers some modification [5]. In contrast, in heavy-ion collisions at typical SIS18 energies, 
i.e., at beam energies around 1 AGeV, the in- medium behavior of vector mesons in compressed nuclear 
matter can be studied, where temperature effects are small and may be neglected. Complementary to 
heavy-ion collisions one can seek for in-medium modifications in reactions of hadronic projectiles jH] or 
real and virtual photons |7| at nuclei, as already at nuclear saturation density sizeable modifications of 
vector mesons are predicted. 

There exists a vastly extended literature on the in-medium modification of hadrons. We mention here 
only the Brown-Rho scaling hypothesis, according to which a mass shift of a vector meson is directly 
interrelated to a change of the chiral condensate jH], and the vector manifestation Pj, the effective 
Lagrangian approach [TU], purely hadronic approaches [H1E2], anc ^ QCD sum rules [T3J EH E3 EH] • 
QCD sum rules [T7] follow the idea of duality (cf. jUI) by relating quantities expressed by partonic 
(quark and gluon) degrees of freedom with hadronic observables. We take here the attitude to assume 
that the partonic quantities are given and examine the QCD sum rule to elucidate the in-medium change 
of the p and uj meson on a common footing. The corresponding current operators, expressed by the 
interpolating quark field operators u and d have the form J? = ^(uj^ l u — d'j f j i d) and = \(u'y fl u + d'y^d), 
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suggesting that the same partonic quantities enter the sum rules. A common treatment of p and to 
mesons is necessary since our schematic transport model studies [HI] show that the relative shift of the 
peak positions and in-medium broadenings must be known to arrive at firm predictions concerning the 
chances to identify both mesons in experiments performed at HADES. 

In a broader context such studies address the phenomenon "mass". Within the standard model, 
masses of quarks and leptons are generated by spontaneous symmetry breaking (in the Higgs mode) in 
the electro-weak sector, while the spontaneous symmetry breaking (in the Goldstone mode) explains the 
features of the hadronic mass spectrum: light Goldstone bosons (71", K, rj with finite masses generated 
by an additional explicit symmetry breaking) and heavier hadronic states emerge j2D|- Within the 
QCD sum rule approach [T7J [TH] the hadron masses can be related to condensates which represent non- 
perturbative features of the QCD vacuum. This is highlighted by the Gell-Mann-Oakes-Renner relation, 
m 2 f 2 = —2m q (qq) (which is actually related to PCAC), and the Ioffe formula, m N = —8Tr 2 (qq)/M 2 (cf. 
21j). Both expressions, which are in leading order, suggest that hadron masses are tightly related to 
the chiral condensate (qq) , which in turn is a measure of the order parameter of the chiral symmetry 
breaking. 

Since (qq) changes with changing temperature j22j and density (cf. Fig. 1 in [22]), one can expect an 
simultaneous change of the hadron masses. This was the very idea of the Brown-Rho scaling hypothesis 
jH] which is modified according to [13]. This promoted very much the physics programme at HADES 
with the motivation to verify, via an observed change of the hadron masses, the chiral condensate's 
change. 

To call a loosely spoken analogy, one is seeking for the QCD analog of the Zeeman/Stark effects: 
an external field (here: strong interaction mediated by surrounding hadrons) changes the excitation 
spectrum of an atom (here: the hadronic excitation spectrum above the QCD ground state). Very often 
the notion "mass shift" is used as short hand notation for a shift of the peak in the spectral function. 

Our paper addresses such mass shifts of p and uo mesons on the basis of the Borel QCD sum rule; 
it is organized as follows: In section II we recapitulate the basics of the QCD sum rule approach. In 
section III we present a combined study of the in-medium behavior of p and u mesons with accounting 
for effects of finite widths. The p — u mass splitting is discussed in section IV. The conclusions can be 
found in section V. 



2 QCD sum rule 

Within QCD sum rules (QSR) the in-medium vector mesons V = p, u are considered as resonances in 
the current- current correlation function 

n) = ij d A x j* x (T jj(ar) jj(0)) n , (1) 

where q^ = (E, q) is the meson four momentum, T denotes the time ordered product of the respective 
meson current operators J^(x), and (■ ■ •)„ stands for the expectation value in medium. In what follows, 
we focus on the ground state of low-density baryon matter approximated by a Fermi gas with nucleon 
density n. We consider isospin symmetric nuclear matter, where the p — u mixing effect is negligible. 

The correlator (JTJ can be reduced to |n^(g 2 ,ra) = Il( v \q 2 ,n) for a vector meson at rest, q = 0, in 
the rest frame of matter. In each of the vector meson channels the corresponding correlator H^ v \q 2 ,n) 
satisfies the twice subtracted dispersion relation, which can be written with Q 2 = — q 2 = —E 2 as 

with nW(0,n) = n( y )(g 2 = 0,n) and IlW(O) = dU ^ q2) \ g ^=o as subtraction constants, and R {v \s) = 
-Imn( y )(s,n)/(7rs). 
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For large values of Q 2 one can evaluate the r.h.s. of eq. (0) by the operator product expansion (OPE) 
leading to 



n( y )(Q 2 



-<* HQ 2 ) + E -m 

i=l 



(3) 



where the coefficients Cj include the Wilson coefficients and the expectation values of the corresponding 
products of the quark and gluon field operators, i.e. condensates. Performing a Borel transformation 
of the dispersion relation (J2J) with appropriate parameter M 2 and taking into account the OPE © one 
gets the basic QSR equation 



°? oo 

nW(0,n) + J dsRW{s) e- s ' M2 = c M 2 + ]T 



(4) 



The advantage of the Borel transformation is (i) the exponential suppression of the high-energy part 
of R v (s), and (ii) the possibility to suppress higher-order contributions to the r.h.s. sum. Choosing 
sufficiently large values of the internal technical parameter M one can truncate the sum in a controlled 
way, in practice at % = 3. 

For the subtraction constants n (y) (0,n) in eq. © we use n (p) (0,n) = n/{AM N ), U^(0,n) = 
9n/(4M7v), which are actually the Thomson limit of the VN scattering processes, but also coincide 
with Landau damping terms elaborated in |2I] for the hadronic spectral function entering the dispersion 
relation without subtractions. For details about the connection of subtraction constants and Landau 
damping term we refer the interested reader to |25j . 

In calculating the density dependence of the condensates entering the coefficients q we employ the 
standard linear density approximation, which is valid for not too large density. This gives for the chiral 
quark condensate (qq) n = (qq)o + > where we assume here isospin symmetry for the light quarks, 
i.e. m q = 5.5 MeV and (qq)o = — (0.24GeV) 3 . The nucleon sigma term is = 45 MeV. The gluon 
condensate is obtained as usual employing the QCD trace anomaly (^G 2 ) n = (^G 2 ) — §Mjy n , where 
a s = 0.38 is the QCD coupling constant and Mjy = 770 MeV is the nucleon mass in the chiral limit. 
The vacuum gluon condensate is (^-G 2 )o = (0.33 GeV) 4 . 

The coefficient C3 in eq. (J3J contains also the mass dimension-6 4-quark condensates (cf. |2S] for 
a recent calculation of corresponding matrix elements) ((?7 A tA a g) 2 ) n , {{u^^X a u){d^X a d)) ni {(qj^\ a q) 
(s7 M A a s)) ra , and ((?7/i7 5 A a g) 2 ) n which are common for p and uj mesons. On this level, p and u mesons 
differ only by the condensate ±2((u7 ^\ a u)(d^^\ a d)) n (cf. [16J), causing the small p — uj mass 
splitting in vacuum [17] . Keeping in mind the important role of the 4-quark condensate [23 EH] for the 
in-medium modifications of the vector mesons, we employ the following parameterization 



<(w 5 A a g) 2 >, 
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ko m q {qq) 
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The parameter k reflects a deviation from the vacuum saturation assumption. (The case k = 1 
corresponds obviously to the exact vacuum saturation JSj-) To control the deviation of the in-medium 
4-quark condensate from the mean-field approximation we introduce the parameter k^. An analog 
procedure applies for the other 4-quark condensates, each with its own k Q and k^, which sum up to 
a parameter k and a parameter k^. As seen in eq. (J3J) and eq. © below, k n parameterizes the 
density dependence of the summed 4-quark condensates; k is adjusted to the vacuum masses. Below 
we vary the poorly constrained parameter kn to estimate the contribution of the 4-quark condensates 
to the QSR with respect to the main trends of the in-medium modification of the vector meson spectral 
function. (Strictly speaking, k and k n differ for p and u mesons due to contributions of the above 
mentioned flavor-mixing condensate; in addition, in medium a twist-4 condensate makes further p and 
uj to differ (T3]. However, the differences can be estimated to be sub-dominant. Therefore, we use in 
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the present work one parameter k^, keeping in mind that it may slightly differ for different light vector 
mesons.) 

Using the above condensates and usual Wilson coefficients one gets as relevant terms for mass 
dimension < 6 and twist < 2 (23 EH] 



f 



c 

Cl 

c 2 
c 3 



f 



7T 



8tt 2 
3m 2 

m q (qq)o + — n + — 



112 ,_ > 2 

— 7T a, k (gg) 
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<^G 2 ) - M« n 
K m g (qq) 
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[ " s : (qq)lQl(2±9) 



1 H ; — — n 



m q (qq} 



(9) 



The terms with A 2 ^ in 02,3 correspond to the derivative condensates from non-scalar operators as a 
consequence of the breaking of Lorentz invariance in the medium. These condensates are proportional 
to the moments of quark and anti-quark distributions inside the nucleon at scale /i 2 = lGeV 2 (see for 
details [13J ) . Our choice of the moments A 2 and is 1.02 and 0.12, respectively. The last line in eq. Q 
stems from the flavor mixing condensate, (u- ■ -ud- ■ -d) n , which has been evaluated with a technique 
similar to ^7]. Qq ~ (9(200 MeV) is a cut-off parameter from momentum integrals. For our purposes 
we can neglect terms related to Qq. 

The value of k in eq. © is related to such a choice of the chiral condensate (qq)o to adjust the 
vacuum vector meson masses. In our QSR we have used kq = 3, obtaining m 



n 



0) = 777 MeV 

close to the nominal vacuum values. The ratio kn/ko m the parameterization (JHJ) is restricted by the 
condition ((?7/iA a g) 2 ) n < 0, so that one gets < < 4 as reasonable numerical limits when considering 
n < n , as dictated by our low-density approximation. 

The case of finite baryon density and temperature has been considered in [2*7j . Here we focus on 
density effects with the reasoning that temperature effects below 100 MeV are negligible. 

To model the hadronic side of the QSR (jlj) we make the standard separation of the vector me- 
son spectral density into resonance part and continuum contribution by means of the threshold 
parameter sy 

,yl l( ' s '" ) e(s v -s) + c e(s-s v ) , (10) 



s,n 



v 



where S^ v '(s,n) stands for the resonance peak in the spectral function; the normalization Fy is unim- 
portant for the following consideration. In vacuum, this ansatz is justified since the time-reversed 
reaction, e + e~ — > V, which is directly related to By\ exhibits a prominent vector meson peak at low 
energies and a smooth continuum at higher energies. The sum rule can be then cast into the form 



/ ds S^ v \s,n) e- fi / M2 



[1 - (l + £) e-/*'] 



c 2 c 3 
M 2 M A 



J ds S^(s,n) s' 1 e~ s / M Cq h 



M- M 4 2M 6 



n^)(o, 

M 2 



n 



Given as such, the sum rule allows to test consistency of a particular model for the spectral function 



as exercised, e.g., in [TQJ . In the zero-width approximation, S^^s^n) oc 5(s 
of eq. (|TT|) becomes simply 



7,2 



m 



n)), the l.h.s. 

In this sense one could also consider the l.h.s. as averaged mass, 
n), keeping in mind that it represents a normalized moment of the hadron strength 
S( V K Then the sum rule eq. (JTTJ determines the parameters m v or thy by the density dependence 



rnyin) 



denoted as rriyi 
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of the condensates, encoded in the coefficients 01,2,3 from eqs. (0 ■ • - EH), and the subtraction constant 
11^(0, n). Without further explication of nothing can be deduced from eq. on in-medium 
mass shifts and broadening (see, however, where, after determining the mass parameter my(n) an 
estimate of the corresponding width is attempted). 



3 Finite width effects 

Our intention is now to evaluate the QCD sum rule eq. (fTTj) by taking into account the finite widths of 
the vector mesons. Note that already in vacuum the p and uj widths differ noticeably, 150.7 MeV and 
8.43 MeV, respectively. This must be reflected in the form of S™. 
Ref. made for the p meson the ansatz 

S(V) ^ n) = ( S - mul)) 2 + (lv(n)T v (s)y (12) 

with the two parameters, my(n) and 7y(n), and some parameterization of the width IV (s). Clearly, 
the one eq. cannot determine the two unknowns my and 7y, rather only the correlation of 7y(my) 
is determined at a given density. In line with the above ansatz ()12|) . one can use a more realistic form 
for the resonance spectral density based on the general structure of the in-medium vector meson 
propagator 

S^\s,n) = - 2 Mn)MS v (s,n) ^ 

(s — m v (n) — Re£y(s, n)) 2 + (yv(n) Im£y(s, n)) 2 

with ReEy(s,n) and Im £y(s,n) as real and imaginary parts of the in-medium vector meson self- 
energy. In the spirit of eq. ()12)) [TE] the meson mass parameter my(n) and the width factor 7y(n) 
become density dependent in nuclear matter to have a degree of freedom for the "request" of the 
QSR. This dependence is determined by the QCD sum rule eq. (fTTj) and mainly governed by the QCD 
condensates. (An analogous approach with 7 P = 1 was used in [Hj.) The in- medium vector meson mass 

2 

is determined by the pole position of the meson propagator, i.e., m v = m v (n) + ReEy(s = m v (n),n). 
Within the linear density approximation the vector meson self energy is given by 

Z v (E,n) = Y^ C {E) - n T VN (E) , (14) 

where E = y/s is the meson energy, Sy c (i?) = Y>y(E,n = 0), and Tvn(E) is the (off-shell) forward 
meson-nucleon scattering amplitude in free space. The renormalized quantity £™ c is summarized in 

2 

the Appendix A in [29|. For the uj meson we absorb as usual ReS™ c in and put simply ImS™ c = 

—m^T^Qi^E — 3m n ) with the vacuum values of mass m w and width r w . 

The described framework is well defined, supposed Ty^ is reliably known. 1 Unfortunately, the 

determination of Ty^ is hampered by uncertainties (cf. results in ^0] arid J2])- ImTyAr is more directly 

accessible, while ReTy^ follows by a dispersion relation with sometimes poorly known subtraction 

coefficients. Since our emphasis here is to include the collision broadening and other finite width effects 

2 

in the spectral function, we absorb, as intermediate step, ReTyjv in m v {n) thus neglecting a possible 
strong energy dependence. In such a way, the uncertainties of ReTyjv become milder since my (n) is then 
mainly determined by the QSR. Neglecting the energy dependence of ReTy^v one discards a potentially 
rich structure of S^ v \ such as, for instance, a double peak structure obtained in jT2j for the uj meson 
or in ^1] for the p meson. Afterwards the importance of ReTy^r is checked. 



1 Then, ~fy = 1 and my(n) = my(0), and the QSR acts merely as consistency check, as mentioned above. 
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We take the needed scattering amplitude Tvn(E) for p and ui mesons from results of the detailed 
analysis of pion- and photon-nucleon scattering data performed recently in ^2] on the footing of the 
Bethe-Salpeter equation approach with four-point meson-baryon contact interactions and a unitary 
condition for the coupled channels. 




m p [GeV] mco [GeV] 

Figure 1: Results of the QCD sum rule evaluations for p (left panels) and to (right panels) 
at n . T VN is discarded. 




0.4 0.5 0.6 0.7 0.8 q.6 0.7 0.8 0.9 1.0 



m p [GeV] m a) [GeV] 

Figure 2: As in Fig. 1 but with ImTyN from \V2\ . 

To test the importance of this particular form of Tvn we proceed to evaluate the sum rule 2 in three 
steps: (i) first neglect Tvn at all, (ii) include only ImTyN, and (iii) include both ImTyN and ReTyN 

2 At a given baryon density n the continuum threshold sy is determined by requiring maximum flatness of 
my(n; M 2 , sy) as a function of M 2 within the Borel window M^ in ■ ■ ■ M 2 nax . The minimum Borel parameter M^ in 
is determined such that the terms of order 0(M~ e ) on the OPE side eq. (J2J contribute not more that 10%. Selecting such 
sufficiently large values of M^ in suppresses higher-order contributions in the OPE eq. 10} and justifies the truncation. 
Typically, M^ in (10%) is in the order of 0.6 GeV 2 . The values for M^ ax are roughly determined by the "50% rule", i.e., 
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Figure 3: As in Fig. 2 but with both ReTyjv and \mTvN from [T2j . 



as well. The results are exhibited in Figs. 1-3. Compared to vacuum, the in-medium changes of the 
r.h.s. of eq. (fTTf require for the p meson more strength at low energy which is accomplished, for the 
given parameterization, by a down-shift of the peak position or a larger width. With increasing 7 the 
masses are shifted to larger values. This can be understood in the following way: A larger width gives 
more contribution to the integrals on the l.h.s. of eq. (fTTj) at lower energy; this is compensated by an 
up-shift of the peak position. Note that 7 = reproduces the zero-width approximation. 

One observes in Figs. 1 - 3 a strong sensitivity on the density dependence of the 4-quark condensate. 
The sensitivity against variations of kn is larger than that of 7. (The displayed range of 7 covers an 
enormous range of the widths.) Remarkable is the tendency of a down-shift of the p mass and a related 
up-shift of the uj mass when restricting the width parameters jv to smaller values. The overall pattern 
seen in figs. 1 - 3 is quite robust. The numerical details, however, change under variations of Tvn- To 
have a fix point let us consider kn = 3 (which will later turn out as relevant value). The peak position 
of the p meson is stable with respect to variations of Tvn- This may be attributed to the large vacuum 
width of p; changes of the amount of T VN by 7 cause only small changes of the p peak position. The 
u meson, in contrast, seems to sit at some borderline: the variations of Tvn can cause an up-shift or a 
down-shift. Inclusion of lmT V N alone results in a tiny effect, while ReT V N pushes the peak further up. 
Larger values of kn are required to get the uj meson's peak position down-shifted. 

In contrast to the universal scaling hypothesis [8. there is no unique in-medium behavior of p and uj 
mesons. Rather, the strikingly different values of the subtraction constants lly(0, n) cause the different 
behavior of p and uj mesons, as stressed in |24j. 3 



the continuum part of the hadronic side must not contribute more than 50% to the total hadronic side to be sufficiently 
sensitive to the resonance part. According to our experience |27l 1281 l2"§| . my is not very sensitive to variations of Af^ ax . 
We can, therefore, fix the maximum Borel parameter by M^ ax = 1.5 (2.4) GeV 2 for the ui (p) meson, in good agreement 
with the " 50% rule" . The sensitivity of the results on these choices of the Borel window is discussed in ■ Flat curves 
mv(n; M 2 , sy) within the Borel window represent a prerequisite for the stability of the sum rule analysis. 

3 Thc detailed analysis of further terms in the OPE side of the sum rule will be relegated to a separate study. 
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1200 




Figure 4: The peak positions of the uj (upper curve) and p meson (lower curve) as a function 
of the parameter kn- Same choice of parameters as in Fig. 2 for 7 = 1. The arrows depict 
the vacuum masses. 



4 p — uj mass splitting 



Inspecting Figs. 1-3 seems to point to a lacking predictive power of the QSR. However, if one is 
interested in the peak position of the spectral function (usually called "in-medium mass"), and not in 
the very details of the shape of S^ v \ and if one assumes that the essential features of the in-medium 
broadening are sufficiently accurate described by ImTyN, one has to put j v = 1. Then one arrives at 
the in-medium mass splitting of p and uj mesons as displayed in Fig. 4. The mass splitting is surprisingly 
large already at normal nuclear density no- The magnitude of the mass splitting can be traced back to 
the particular form of ImTyN from ^2]- (Using ImTyN from PU| results in a smaller splitting.) The 
driving force of the p - uj mass splitting, however, is the difference of n w (0,ri) and IP(0,n). 

The p and uj mass shifts are determined by the still unconstraint parameter kn- Once one of the 
vector mesons's peak position is experimentally determined the other one is fixed by the correlation 
displayed in Fig. 4, up to some uncertainty caused by the actual lmT VN or the freedom in j v - 

A presently running experiment [?J measures the reaction 7 + A — > X + uj with identifying the uj via 
the 37 Dalitz decay. First data analysis [7j seems to exclude any up-shift, and even a weak down-shift 
of uj strength is compatible with the data. If this result gets confirmed it would require, within the 
present approach, that k n > 3, i.e., a strong density dependence of the 4-quark condensate. In other 
words, this would give a first direct evidence for a change of a condensate. Even a null effect for uj 
mesons decaying inside the target nucleus would require a strong density dependence of the 4-quark 
condensate, as evidenced by Fig. 4. In fact, at n and for = 3 the 4-quark condensate drops to 60% 
of its vacuum value. 
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5 Conclusion 



In summary, we present here an analysis of the QCD sum rule for the in-medium behavior of p and u 
mesons. Truncating the "QCD side" of the sum rule at mass dimension 6 and twist 2 we find a strong 
sensitivity of the p and uj peak positions on the density dependence of the 4-quark condensate. Assuming 
that we have at our disposal a sufficiently realistic description of the collision broadening, described in 
low-density approximation by the off-shell forward meson- nucleon scattering amplitude ImTyN, the in- 
medium mass shifts of p and uj mesons are related essentially to one parameter. Experiments identifying 
the uj decay in nuclear matter then constrain this parameter, thus fixing also the in-medium p mass, up 
to uncertainties inherent in the sum rule approach. If an up-shift of uj strength can be experimentally 
excluded, the present QCD sum rule analysis points to a strong reduction of the 4-quark condensate 
already at nuclear saturation density. 

On a quantitative level there is some uncertainty caused by the actual form of IhiTVat, in particular 
for the p meson which influences the uj — p mass splitting. Poorly known twist-4 condensates modify 
further this splitting and deserve additional investigations. Given the importance of the 4-quark con- 
densate, one could be afraid on the influence of higher order condensates. Here, the hope is that an 
appropriate Borel window suppresses these higher orders. We also found some quantitative changes 
when including explicitly ReT VN in the spectral function. But the overall pattern of the in-medium 
modifications of p and uj mesons is stable. 

In contrast to p and uj mesons, which are insensitive against changes of the genuine chiral condensate, 
(qq), the meson depends sensitively on the chiral strange condensate, (ss), and only very weakly on 
the 4-quark condensate. 

With respect to the future accelerator project SIS200/300 at GSI an extension of the present ap- 
proach to the in-medium behavior of D mesons is challenging. 

One of the authors (B.K.) thanks Amand Faessler for inviting him to the International School of Nuclear 
Physics, Erice (It), Sep. 16 - 24, 2003, where the material of this article has been presented. Many useful 
discussions with S. Leupold, M. Lutz, V. Metag, U. Mosel, and Gy. Wolf are gratefully acknowledged. The 
work is supported by BMBF 06DR121 and GSI. 
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